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Low temperature specific heat and thermal conductivity measurements on the ambient pressure 
heavy fermion superconductors Celrlus and CeCoIns reveal power law temperature dependences 
of these quantities below Tc. The low temperature specific heat in both Celrlns and CeCoIus 
includes terms, consistent with the presence of nodes in the superconducting energy gap. The 
thermal conductivity data present a T-linear term consistent with the universal limit (Celrlns), and 
a low temperature T"^ variation in the clean limit (CeCoIns), also in accord with prediction for an 
unconventional superconductor with lines of nodes. 
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Unconventional superconductivity has been an ac- 
tive area of research for several decades, ever since 
the discovery of the first heavy fermion superconductor, 
CeCuaSia fill. The presence of strong magnetic interac- 
tion between 4f moments and itinerant electrons in this 
class of compounds allows the possibility of non-phonon 
mediated coupling between superconducting quasipar- 
ticles, a signature of unconventional superconductivity, 
and a superconducting order parameter with lower sym- 
metry than that of the underlying crystal lattice. Soon 
after the discovery of superconductivity in CeCu2Si2, sev- 
eral uranium-based heavy fermion superconductors were 
discovered. The presence of a double transition in UPta 
immediately identified this compound as an unconven- 
tional superconductor ||]. Subsequent observations of 
power law temperature dependences of the specific heat 
and thermal conductivity have been instrumental in iden- 
tifying UPta (and other U-based heavy fermion supercon- 
ductors) as unconventional. 

Until very recently, among Ce-based heavy fermion 
compounds only CeCu2Si2 was shown to superconduct 
at ambient pressure. Other ThCr2Si2-based compounds 
require application of significant pressure (on the order of 
20 kbars) before they exhibit superconductivity. These 
include CeCu2Ge2 [|, CePd2Si2 |,|, and CeRh2Si2 [|. 
Recently, cubic Celna was also shown to superconduct 
under pressure of about 25 kbar |Q, with superconduc- 
tivity mediated by magnetic interactions |^ . 

Quantitative thermodynamic and heat transport mea- 
surements, that probe the nature of the superconducting 
gap, are a lot simpler at ambient pressure. CeCu2Si2 
seems to be a good candidate for such investigations. 
However, it is a rather complicated system, in which 
very small changes in stoichiometry can change the 
ground state from superconducting to antiferromagnetic. 
Recently discovered ambient pressure superconducting 
Celrlns H and CeCoIns H do not display such complica- 



tions. Therefore, they present an uncommon opportunity 
to study the superconducting ground state of Ce-based 
heavy fermion compounds. In this Letter we present the 
results of low temperature specific heat and thermal con- 
ductivity measurements on both Celrlns and CeCoIns 
which show power law temperature dependences and 
present a strong proof of the unconventional nature of 
the superconductivity in these compounds. 

The details of sample growth and characterization are 
described in Refs. Large plate-like single crystals, 

up to 1 cm long, are grown from an excess In flux. Their 
quasi-2D tetragonal structure can be viewed as layers of 
Celna separated by layers of Irln2 (or CoIn2). Therefore, 
we can treat Celua as the parent compound for Celrlns 
and CeCoIus. Both Celrlns and CeCoIus are ambient 
pressure superconductors, with transition temperature 
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FIG. 1. Specific heat of Celrlns, as (C - Csch)/T vs. T. 
(□) H = kC, (o) H = 5 kC. Dashed line is 7 = C/T = 0.68 
J/molK^; solid line is a fit to H = kC data for 85 mK <T < 
0.2 K. Inset; C(T), (□) H = kC, (o) H = 5 kC. Solid line is 
a sum of 7T and nuclear Schottky Csch for H = 5 kC data. 
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FIG. 2. (C - Csch)/T of CeCoIns at H = kG. Solid line 
is a fit for 0.1 K < T < 0.7 K. Inset: Specific heat of CeCoIng: 
Solid line is In nuclear quadrupolar Schottky anomaly Csch- 

Tc of 400 mK and 2.3 K and specific heat jump at su- 
perconducting phase transition AC/jTc of 0.76 and 4.3 
respectively. 

Fig. |l| shows specific heat of two samples of Celrlus in 
zero field and in magnetic field of 5 kG (H || c), which 
suppresses the superconducting state to T = K. The 
solid line in the inset represents the sum of the T-linear 
term expected for a well developed heavy fermion state 
and the fn nuclear quadrupolar Schottky anomaly Csch, 
and it agrees well with the 5 kG normal state data for 
Celrlus . The five doubly degenerate energy levels used in 
calculating Csch (In nuclear spin / = 9/2) were measured 
directly via NQR experiment . 

The main body of Fig. |^ shows specific heat divided 
by temperature after C'sch was subtracted. The solid line 
in the main body of Fig. ^ is a linear fit to the data in 
the superconducting state of Celrlus below 0.2 K {Tc/2). 
The vertical intercept gives the coefficient of the T-linear 
term 70 = .110 ± .010 J/molK^ (a second sample shown 
in the inset gives 70 = .010 ± .008 J/molK^ between 
200 mK and 85 mK). Low temperature T-linear contri- 
bution in specific heat of unconventional superconduc- 
tors is commonly attributed to the impurity band that 
forms in the linear node(s) of the superconducting en- 
ergy gap. The sample with larger 70 has lower Tc of 
0.38 K, compared to 0.4 K for the sample with lower 
7o, qualitatively consistent with the impurity band ori- 
gin of the linear-T term in specific heat. The term 
in electronic specific heat in Celrlns is an indication of 
the presence of lines of nodes in the energy gap, as in a 
number of other unconventional superconductors, includ- 
ing both cuprates YBCO |ll| and LSCO |l|] and heavy 
fermion UPd2Al3 ||l^]. The coefficient of the T-squared 
term in the specific heat, a = 3.56 ± 0.06 J/molK"^, of 
Celrlus is several orders of magnitude greater than that 
of the HTS, reflecting the heavy- fermion nature of this 



compound. 

Specific heat of CeCoIns in zero field is shown in Fig. ||. 
The inset of Fig. |2| shows the data up to 3 K together 
with calculations of the In quadrupolar nuclear Schottky 
anomaly (solid curve) based on the energy levels deter- 
mined via NQR studies of CeCoIns Q . The good agree- 
ment between calculation and data leaves little doubt 
that the low temperature anomaly is due to In nuclear 
quadrupolar moments. In the main body of Fig. | we 
show specific heat divided by temperature after subtract- 
ing the In nuclear Schottky contribution. The zero field 
data are well described by C/T = 0.04 J/molK^ + 0.25r 
J/molK'^ for 95 niK < T < 400 mK suggesting the pres- 
ence of impurity band and line(s) of nodes in the energy 
gap, as in the case of Celrlus. The low temperature 
upturn in the data may be related to the low tempera- 
ture entropy revealed when 5 T field is applied to sup- 
press superconductivity in this compound. Crystallo- 
graphically, Celrlns and CeCoIus are closely related com- 
pounds. The similar temperature dependences of their 
specific heats points (gratifyingly) to the possibility of 
the same symmetry of their superconducting order pa- 
rameters. 

We now turn attention to thermal transport in Celrlns 
and CeCoIus. Thermal conductivity k of Celrlus divided 
by temperature is displayed in Fig. ^ for 2 K < T < 33 
mK. k/T reaches a maximum at Tc = 0.4 K, begins to 
fall below Tc, and becomes linear in temperature below 
Tc/2 = 200 mK. The solid line is a fit k/T = .46 ± .07 
W/K^m +(4.10 ± .06)T W/K^m for T < 200 mK. The 
dashed line in Fig. ^ is an upper limit estimate of the 
phonon thermal conductivity (with boundary scattering 
only) K ^ i/3T3(u)Ao. Here PT^ = 15.5T^ is 
the phonon specific heat of the non-magnetic analogue 
Lalrln5 fisll , Aq is the phonon mean free path, and (v) is 
the mean phonon velocity [^,0. On the basis of this 
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FIG. 3. Thermal conductivity of Celrlns. Solid line is a 
linear fit to the data for T < 0.2 K = Tc/2. Dotted line is an 
upper limit estimate for the phonon thermal conductivity. 
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estimate we can conclude that the heat transport below 
Tc is dominated by quasiparticles. 

The presence of a T-linear term in k is expected for an 
unconventional superconductor with line(s) of nodes in 
the energy gap on the Fermi surface. It arises due to the 
appearance of the impurity band at these nodes. At first 
sight the magnitude (0.46 W/K^m) of the coefficient of 
the linear term seems to be inconsistent with the results 
of the specific heat measurements, which indicate a very 
clean sample with low impurity concentration. However, 
Graf et al. p8| has shown that the low temperature limit 
of thermal conductivity in unconventional superconduc- 
tors with a variety of order parameter symmetries that 
have lines of nodes is universal (independent of the con- 
centration of scattering impurities), 

— = — NfVf X — . (1) 

Here Nj is the normal density of states, is a Fermi 
velocity, Aq is a maximum of the energy gap at T — 0, 
fi is the slope parameter which describes the rate of the 
increase of the energy gap away from the node on the 
Fermi surface, and a is 1 for 3D order parameter with 
lines of nodes and - for the 2D case (as in the d-wave 
state of the high- temperature superconductors (HTS)). 
In the rest of this Letter and in variety of calculations of 
Ref. /i was taken to be equal to 2. 

We obtain the normal state density of states Nf, Fermi 
velocity v / , and mean free path Itr (or scattering life time 
r) at Tc following the BCS-based analysis of Refs. ||l|,|o|. 
The input parameters for this analysis are the linear coef- 
ficient of the specific heat 7, the derivative of the critical 
magnetic field at T^, and the normal state resistivity 
at Tc po (21). We obtain Nf = 1.13 x 10^6 (ergcm^)"!, 
Vf = 7.65 X 10^ m/s, effective mass ?7ie// ~ 140me, and 
Itr = 1350 A at Tc. These numbers are to be com- 
pared with those obtained from dHvA experiments [^ , 
which have observed We// of up to 45me, Fermi veloc- 
ities down to 6.3 x 10^ ■j-, and quasiparticle mean free 
paths between 58 A (for the lightest band) and 4500 A 
(for the heaviest band observed). Overall, the results of 
the analysis above compare favorably to those of dHvA 
experiments. We also obtain zero temperature coherence 

length ^0 = \/ 27THatlT=6 ) ^ ^41 A from the H-T phase 
diagram and the ratio ltr/£,a = 5.43, which shows 
that Celrln5 is indeed in the clean limit. The energy gap 
is then Ao/ks = 0.746 K, and Ao/ksTc = 1.865, close 
to BCS value of 1.763. 

We can now estimate the universal limit via Eq. ^, 
obtaining y |j,^q = 1-06 W/K^m for a 3D order param- 
eter. If in addition we take Ag = 2.1AkBTc, a weak cou- 
pling value for the d-wave order parameter in the clean 
limit, we obtain yj^^Q = 1.2 W/K^m for a 2D d-wave 
order parameter. Given the fact that the Orlando anal- 
ysis assumes a spherical Fermi surface, while LDA 



band structure calculations |2^] reveal a very compli- 
cated Fermi surface with contributions from three differ- 
ent bands, an uncertainty in estimated of about 
a factor of two seems reasonable. Therefore the large ob- 
served yjy^Q — 0.46 W/K^m can be explained on the 
basis of the theory of Graf et al. |l^, and is additional 
strong proof of the presence of line(s) of nodes in the gap 
of Celrln5. 

Finally, Fig. ^ shows thermal conductivity of CeCoIns 
on a log-log plot. The data are striking in several ways. 
First, upon the sample entering the superconducting 
state, thermal conductivity displays a sharp kink and 
rises from 2 W/Km at Tc — 2.3 K to a maximum value 
of 5 W/Km at 0.71 K (T/Tc = 0.3), an increase of a fac- 
tor of 2.5. It is more instructive to consider the value of 
k/T, which is proportional to the product of the density 
of normal quasiparticles and their life time. k/T grows 
from 0.86 W/K^m at Tc ^ 2.3 K and reaches the value 
8.7 W/K^m at T = 0.46 K (T/Tc = 0.2). In spite of the 
strong reduction in the number of normal quasiparticles 
expected at T = 0.2Tc, k/T actually grows by more than 
an order of magnitude! An increase of k/T was also ob- 
served in several HTS, where the separation of electronic 
and phonon contribution, and therefore identification of 
the origin of the peak in k, is rather complicated p^ , p5[ . 

In the case of CeCoIns the situation is much more 
straightforward. The dotted line in Fig. ^ is an upper 
limit estimate of the phonon thermal conductivity Kph in 
CeCoIns. At the temperature of the maximum in k the 
phonon upper bound is a factor of 15 below the measured 
value. The dashed line for T > Tc is a Widemann- Franz 
estimate of electronic thermal conductivity in the normal 
state Kg = LqT/p, with the normal state resistivity p 
from Ref. ^ and the Lawrence number Lq = 2.44 x lO"'^ 
Wr2K~^. In the normal state the thermal conductivity 
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FIG. 4. Thermal conductivity of CeCoIns. Solid line is a 
power law fit for T < 100 mK < Tc/20. Dotted line is an 
upper limit estimate of Kph- Dashed line is Ke for T > T^. 
Inset; k/T vs. for T < 0.1 K. 
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of CeCoIn5 is dominated by electrons, which contribute 
over 90% of the experimental k at Tc- We can thus con- 
clude that all of the peak below Tc in CeCoIns is due to 
electrons, and may be viewed alternatively as supporting 
the mostly electronic origin of the peak in k in HTS as 
well. This peak may bear a direct relation to the ap- 
pearance of the low temperature rise in the normal state 
specific heat at 7? = 5 T, which compensates the entropy 
from a very large specific heat jump ^P- — 4.3 at Tc in 
zero field This large ^ suggests that the fluctua- 
tions in the normal state, responsible for the low temper- 
ature tail in C in the 5 T data, are redistributed up to 
Tc when CeCoIns enters the superconducting state. Such 
fluctuations also scatter electrons, and their suppression 
in the superconducting state would increase k/T. 

We estimated the universal limit of k/T in CeCoIns 
via an analysis similar to that for Celrlns ||2^] . The dirty 
limit term contributes only 6% to H'^2i ^tnd Itr/^o ~ 14 
at Tc, clearly indicating that CeCoIns is also a clean su- 
perconductor . Using the values of Aq and a = 1 (for 
a 3D superconductor with lines of nodes) in Eq. |^, we 
obtain for CeCoIng k/T\^^„ ^ 0.19 W/K^m. With the 
d-wave energy gap Aq/Zcb — 2.14Tc and a = 4/7r for a 
2D superconductor, we obtain K,/T\rp^Q k, 0.1 W/K^m. 
The inset of Fig. | shows k/T vs. below 100 mK. The 
minimum value measured is 0.2 W/K^m, which suggests 
that the sample is close to universal limit. Future mea- 
surements to investigate universal limit of CeCoIn5 are 
planned. 

For temperatures below 100 niK (T/Tc « .043) and 
down to the lowest temperature measured of 33 mK (« 
1.5% of Tc) thermal conductivity follows a power law be- 
havior. The straight line in Fig. ^ is a fit K = ar3-37 to 
the data. This is close to T^ behavior predicted for un- 
conventional superconductor in a clean limit with line(s) 
of nodes |l^. Therefore the heat transport in CeCoIns 
between 33 mK and 100 mK can be described as that 
of an unconventional superconductor in the clean limit, 
with an impurity band width of less than 30 mK. An 
upper limit estimate of the impurity concentration that 
gives rise to such narrow impurity band Q is rump = 20 
ppm in the unitary scattering limit and nimp = 20 ppt 
in the Born limit. The true impurity concentration is 
probably between these limits, and closer to the unitary 
scattering limit [||. 

In summary, both Celrlns and CeCoIns display power 
law behavior in specific heat and thermal conductiv- 
ity, indicative of unconventional superconductivity. T- 
squared terms in specific heat imply existence of lines 
of nodes in both superconducting energy gaps. Thermal 
conductivity also supports unconventional order param- 
eters in these compounds, which manifests itself in (1) 
large k/TI^^q = 0.46W/K^m in Celrlns consistent with 
the universal limit estimates, (2) large increase in elec- 
tronic thermal conductivity below Tc = 2.ZK in CeCoIns 



reminiscent of the behavior of HTS, (3) close to T^ behav- 
ior between 33 mK and « IQQmK in CeCoIn5, predicted 
for unconventional superconductors in clean limit. 
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